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Emerging evidence suggests that our memories for recent events depend on a dynamic interplay between 
multiple cortical brain regions, although previous research has also emphasized a primary role for the 
hippocampus in episodic memory. One challenge in determining the relative importance of interactions 
between multiple brain regions versus a specific brain region is a lack of analytic approaches to address this 
issue. Participants underwent neuroimaging while retrieving the spatial and temporal details of a recently 
experienced virtual reality environment; we then employed graph theory to analyze functional connectivity 
patterns across multiple lobes. Dense, large-scale increases in connectivity during successful memory 
retrieval typified network topology, with individual participant performance correlating positively with 
overall network density. Within this dense network, the hippocampus, prefrontal cortex, precuneus, and 
visual cortex served as "hubs" of high connectivity. Spatial and temporal retrieval were characterized by 
distinct but overlapping "subnetworks" with higher connectivity within posterior and anterior brain areas, 
respectively. Together, these findings provide new insight into the neural basis of episodic memory, 
suggesting that the interactions of multiple hubs characterize successful memory retrieval. Furthermore, 
distinct subnetworks represent components of spatial versus temporal retrieval, with the hippocampus 
acting as a hub integrating information between these two subnetworks. 



Spatial and temporal details constitute critical components of our memories for recently experienced events, 
termed episodic memories. Both knowledge about where we were and approximately when things hap- 
pened during our day are important cues for remembering what happened to us at an earlier event. Current 
literature strongly supports the involvement of the hippocampus generally in storage and retrieval of episodic 
memories', specifically in coding both spatial and temporal elements of recent events^ Yet, past studies have also 
implicated regions outside of the medial temporal lobes in episodic memory, such as the prefrontal cortex, and the 
relative importance of this and other cortical areas to episodic memory remains under debate''**. In support of the 
participation of neocortical regions to episodic memory, lesion studies implicate the prefrontal cortex and 
posterior parietal regions in episodic memory functions' '". In particular, lesions to prefrontal cortex have also 
been linked with more specific deficits in temporal order processing", while lesions to parietal/retrosplenial areas 
are linked to more specific deficits in spatial processing'-. Thus, the relative contributions of cortical areas outside 
of the hippocampus to episodic memory, and the exact role of these areas - including the hippocampus - in 
coordinating spatiotemporal information from disparate neocortical areas, remains unresolved. 

Examining successful compared to unsuccessful retrieval of spatiotemporal details, we tested three hypotheses 
regarding large-scale brain networks using functional magnetic resonance imaging (fMRI). Correct retrieval 
might involve greater connectivity to a primary region (e.g., the hippocampus), suggesting it depends principally 
on a single brain area. Alternatively, correct memory retrieval could involve greater increases overall in functional 
connectivity compared to incorrect retrieval, with little contribution from specific regions. Finally, correct 
memory retrieval might involve some combination of the two (overall increase in connectivity and overall greater 
connectivity to multiple areas). This final hypothesis is most consistent with evidence from previous research 
because it incorporates and emphasizes the roles of additional areas (such as prefrontal and parietal cortices) in 
episodic memory retrieval. Our experimental paradigm involved participants retrieving both spatial and tem- 
poral order information independently from a recently acquired experience, thus enabling us to contrast these 
components of episodic memory. To further elucidate these neural correlates, based on past literature^ '-' ''', we 
propose that retrieval of these separate contextual details could involve distinct subnetworks within a larger 
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episodic memory retrieval network, involving a limited but distinct 
set of overlapping hubs (i.e., the hippocampus). Or, spatial and tem- 
poral order retrieval might involve largely overlapping networks, in 
which case we would not find differences when contrasting spatial 
versus temporal order retrieval. To address these competing hypo- 
theses, we chose seed brain regions from 13 bilateral brain areas 
suggested in past work to be involved in some form in episodic 
memory literature but also those related to other processes, such 
as attention and working memory. These areas formed the basis for a 
graph theory analysis in which we directly compared both correct 
and incorrect memory retrieval and spatial and temporal retrieval 
more generally. 

Results 

Participants first navigated a virtual reality environment in which 
they learned both the spatial layout of the city and the temporal order 
of deliveries of a passenger to different stores (Fig. lA). Following 
encoding, we probed participants' contextual memory for their vir- 
tual reality experience (Fig. IB) while undergoing fMRI. Questions 
consisted of a reference store followed by two probe stores; partici- 
pants indicated which of the two probe stores was closer in either 
spatial or temporal distance from the reference store. Behavioral 



accuracy was 64.4 ± 9.7% for spatial and 61.5 ± 9.7% (mean ± 
s.d.) for temporal questions. Importantly, we found no difference 
in performance between conditions (paired t test, t(15) = 0.954, p 
= 0.355; for a complete description of the behavioral design and 
implementation, see Ekstrom et al., 20 11). Analyzing the fMRI data, 
we first focused on two fundamental contrasts using functional con- 
nectivity and graph theory: 1) correct versus incorrect retrieval (i.e., 
trials on which the participants correctly judged the spatial or tem- 
poral distance compared to incorrect trials) and 2) spatial compared 
to temporal retrieval trials. 

Successful versus unsuccessful memory retrieval. To create our 
retrieval networks, we used a beta time series analysis", which 
assumes two regions are coupled during a specific condition of a 
task if the activity of both is significantly correlated across trials. In 
other words, two regions are functionally connected if their 
responses fluctuate similarly under the same condition. To 
characterize differential connectivity between conditions (correct 
versus incorrect contextual retrieval), we took the set difference 
between two networks, i.e., connections that were exclusive to 
either correct or incorrect contextual retrieval (labeled here as 
correct - incorrect and incorrect - correct, see Methods for 




B 



Temporal Context Question 
"Which of the following stores is closer in 
delivery order to the store you just saw?" 



Item Recognition Question 
"Have you seen the following 
store before?" 




Spatial Context Question 
'Which of the following stores is nearer in 
distance to the store you just saw?" 



Figure 1 | Description of Experimental Design. (A) Participants (n = 16) navigated a virtual city where they picked up a passenger (depicted in the left 
panel) and delivered the person to a series of stores; they encoded the spatial layout of the city and the temporal order of deliveries. The middle panel 
provides an example of a store the participant visited. The right panel shows an aerial view of the virtual environment with the stores distributed unevenly 
around a circle. (B) After encoding, participants underwent functional imaging whUe performing the spatiotemporal retrieval task shown here. 
Participants were first presented with a single reference store and asked to indicate whether they saw this store before. After responding, they were 
subsequently presented with the spatiotemporal contextual retrieval question. Participants were shown two additional store pictures and asked to make a 
distance judgment regarding the spatial or temporal proximity to the reference store. 



SCIENTIFIC REPORTS | 4:6431 | DOI: 1 0. 1 038/srep0643 1 



2 



analytic and statistical details). When participants accurately judged 
the spatial or temporal distance of probe questions compared to 
incorrect trials, we found a significant overall increase in large- 
scale network connectivity, indexed by the number of connections 
(also termed edges) between our regions of interest (also termed 
nodes) (x^ (1) = 63.75, p < 0.0001, Fig. 2A,B,C). This result 
suggests that the increase in coordinated activity of multiple 
distributed regions underlies successful retrieval of spatiotemporal 
contextual details. 

Correct retrieval was not only indexed by the density of network 
connections but also by the precise organization of its edges to par- 
ticular nodes. Specifically, we investigated regional importance using 
a combination of two metrics: node degree and betweenness cent- 
rality. Node degree captures the local connectivity of a node, while 
betweenness centrality describes more global aspects of the network 
by tracking the tendency for the connections from other nodes to 
travel through a given node. Previous studies have found these two 
graph theory metrics to be reliable predictors of functional hubs"*"^". 
We found that several areas within the medial temporal lobes (MTL), 
frontal and parietal lobes, and the visual cortex, exhibited signifi- 
cantly greater levels of betweenness centrality than other nodes in 
the network for correct compared to incorrect retrieval trials. 
Figure 2D shows graphs for brain hubs displaying high betweenness 
centrality, with chi-squared tests indicating significant differences 
between the two conditions. Tests were Bonferroni corrected at p 
< 0.05 to account for multiple comparisons between regions (all 
Pcorrected < 0.05; Calcaruie sulcus (Calc): (1) = 93.86; anterior 
hippocampus (aHPC): (1) = 45.92; posterior hippocampus 
(pHPC): (1) = 24.66; inferior frontal gyrus, pars opercularis 
(IFOp): (1) = 15.38; inferior frontal gyrus, pars triangularis 
(IFTr): (1) = 11.00). Two other hubs, precuneus and parahippo- 
campal gyrus, showed higher betweenness centrality for correct - 
incorrect but at an uncorrected threshold of p < 0.05. Thus, our 
analyses supported the idea that correct contextual retrieval involved 
greater overall connectivity across multiple cortical areas, but that 
this connectivity was most marked at specific brain hubs, supporting 
our third hypothesis. 

To rank the relative importance of each hub to correct spatiotem- 
poral retrieval, we organized each node according to its total degree 
(i.e., total number of connections, see Methods; for a comparison of 
betweenness centrality and node degree, see Fig. S3). Consistent with 
our earlier analysis of betweenness centrality, numerous different 
edges were involved in correct spatiotemporal retrieval, with hip- 
pocampus, visual cortex, and inferior frontal gyrus at the top of the 
hierarchy (Fig. 2 A, right). Overall, the total connectivity in the MTL 
(aHPC, pHPC, and pPHG) was greater than that in any other brain 
region (MTL: 51; next highest, visual cortex: 31, %^ (1) = 4.88, p = 
0.027), suggesting a prominent role of this area in successful retrieval 
of spatiotemporal contextual details (Fig. 2A). Interestingly, while 
successful retrieval resulted in numerous hubs and many more func- 
tional connections than incorrect retrieval, one hub did characterize 
the incorrect memory network: inferior frontal gyrus, pars orbitalis 
(IFOr). As shown in the node hierarchy in Figure 2B and in the bar 
graph in Figure 2D, this region of prefrontal cortex had increased 
participation during incorrect rather than correct retrieval {y^ iX) = 
20.02, Pcorrected < 0.05). Past rcscarch has implicated the inferior 
frontal gyrus in suppression of prepotent responses', and it is intri- 
guing to consider that heightened influences of this subarea during 
retrieval could indicate inhibition of other circuits relevant to correct 
retrieval. 

Network density and behavioral performance. Next, we addressed 
whether 1) overall network connectivity, or density, was predictive of 
behavioral performance and 2) if so, did specific hubs contribute 
greater variance than others in terms of explaining this effect. 
Critically, demonstrating a link between memory performance and 



overall network connectivity would provide strong evidence for the 
importance of increases in connectivity in individual participants to 
memory performance. We first employed simple linear regression to 
test if the density of individual networks could significantly predict 
behavioral accuracy. We did this by regressing each participant's 
performance against the total percent connectivity of the 
individual network (see Methods). We found a significant linear 
relationship between participant performance and the overall 
connectivity within their network, which is depicted in the 
scatterplot in Figure 3. These results indicate that the density of 
connections explained 33% of the variance (R^ = 0.33, F(l,15) = 
6.90, p = 0.020, P = 0.22) in participant performance. This 
relationship was not driven by outliers as it was still significant 
using robust regression (p = 0.033). 

We subsequently considered how each individual node contribu- 
ted to this significant linear relationship using a multivariate stepwise 
regression approach. Based on hubs that first showed a significant 
simple regression with performance, a stepwise regression using 
these brain areas (indicated with asterisks in Fig. 3) resulted in a final 
more parsimonious model with MFC, precuneus, hippocampus, and 
IFOr (shown in red) contributing significantly to the model ((3 = 
4.05, 6.83, 6.76, 5.08, -8.83, p = 0.0001, R' = 0.90). These findings 
argue that while overall network connectivity was a significant pre- 
dictor of memory performance, specific hubs drove this effect most 
strongly. Consistent with our earlier analysis, this included connec- 
tivity to the MFC, hippocampus, precuneus (significant for uncor- 
rected values), and IFOr as significant predictors of performance. 
Coefficients for MFC, hippocampus, and precuneus were positive, 
indicating a significant positive linear relationship. IFOr, in contrast, 
displayed a negative coefficient, indicating connectivity at this par- 
ticular node was a negative predictor of performance. These findings 
largely mirror our earlier results, suggesting the importance of hip- 
pocampus, prefrontal cortex, and parts of parietal cortex to success- 
ful memory retrieval. They also provide an important extension by 
showing a direct relationship between connectivity at specific brain 
hubs and individual participant memory performance. 

While our results exhibited a clear relationship between node 
connectivity and behavior, it is reasonable to speculate that the aver- 
age level of activation of each node could also contribute to, or even 
confound, this relationship. For example, nodes that have high levels 
of activation could potentially influence the correlation of the BOLD 
parameter estimates, thus potentially explaining their significant 
functional interactions in the network. To address this idea, we cor- 
related node degree and average parameter estimates (the average of 
parameter estimates across participants for each ROI) and found no 
significant correlation (p = —0.01, p = 0.9614). In addition, adding 
the average parameter estimate for each node to the stepwise regres- 
sion between node degree and behavioral performance contributed 
no significant variance to the overall model. These analyses suggest 
that univariate effects alone are unlikely to be driving our results. 

Differential spatial and temporal connectivity. Comparing 
networks during spatial and temporal order retrieval generally, we 
found striking differences in their topology, as revealed by the two- 
dimensional anatomical plots in Figure 4A,B. During spatial 
retrieval, a greater number of edges were concentrated in posterior 
areas of the brain while during temporal retrieval, connectivity was 
higher in anterior parts of the brain (Fisher's exact test, p < 0.0001, 
Fig. 4C, middle). This was also evident from the adjacency matrices 
(Fig. 4C, left), where the blue and green boxes emphasize fronto- 
temporal connections, and the red box shows connections within 
the entire posterior portion of the brain. 

As shown in Figure 4D, betweenness centrality indicated that the 
SPL and visual cortex had high betweenness centrality for the differ- 
ence between spatial and temporal (spatial - temporal, see Methods) 
compared to the temporal - spatial retrieval networks (all Pcorrected < 
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Figure 2 | Increases in Network Connectivity during Successful Retrieval of Spatiotemporal Details. (A and B) Each ROI (colored circle) is plotted with 
significant connections (pink lines) between the nodes; network during correct context responses shown in A and during incorrect context responses in B. 
Both size and color indicate node degree; a larger radius and warmer colors represent a higher node degree, while a smaller radius and cooler colors 
represent a lower node degree. The ROIs are spatially distributed according to their x and y coordinates in standardized brain space but have been 
collapsed in the z coordinate direction. The right panel displays node hierarchy where anterior nodes (red) and posterior nodes (blue) are plotted in 
ascending order of node degree, with the highest node degree level at the top; nodes belonging to the same horizontal level have the same node degree. The 
abbreviations for the ROIs are as follows: SFG (Superior Frontal Gyrus); MFG (Middle Frontal Gyrus); IFOr (Inferior Frontal Gyrus, Orbitalis); IFTr 
(Inferior Frontal Gyrus, Triangularis); IFOp (Inferior Frontal Gyrus, Opercularis); PrCG (Precentral Gyrus); aHPC (Anterior Hippocampus); pHPC 
(Posterior Hippocampus); pPHG (Posterior Parahippocampal Gyrus); IPL (Inferior Parietal Lobule); SPL (Superior Parietal Lobule); PCN (Precuneus); 
Calc (Calcarine sulcus). (C) Left: The adjacency matrix (right of diagonal: correct - incorrect, left of diagonal: incorrect - correct) depicts the data to 
optimize viewing the intra-lobular connections. The colored boxes indicate intra-lobular connections in the frontal lobe and temporal lobe (blue, green, 
respectively). For successful retrieval, numerous connections exist within the MTL but many are found throughout the neocortex, compared to low levels 
of functional connectivity throughout the brain during unsuccessful retrieval. Right: Comparison of the total number of edges (significant functional 
connections) between conditions, with the correct - incorrect network containing significantly more edges. *p < 0.0001. (D) Bar graphs showing 
betweenness centrality for each ROI (collapsed across hemispheres). Asterisks indicate Pcorrected < 0.05 significant differences using a chi-squared test on 
node connectivity between conditions (i.e. correct - incorrect and incorrect - correct). 
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Figure 3 | Individual Patterns of Network Connectivity Predict Memory Performance. The scatter plot shows a significant linear relationship 
between total connectivity and individual participant behavioral performance (p < 0.05). To determine which ROIs were driving this correlation, each 
node's contribution to overall network density was regressed against individual performance using a multiple linear regression. The initial model in the 
stepwise regression consisted of the nodes labeled with asterisks, and a final more parsimonious model resulted in the nodes highlighted in red (MFG, 
HPC, PCN, and IFG). The beta coefficients from the linear regression and stepwise regression are presented on a line next to the independent variables. 
For example, the network density resulted in a coefficient of 0.20 for the linear regression. 



0.05, SPL: (1) = 15.41 and primary visual cortex: (1) = 20.15). 
Conversely, MFG and IFG nodes, located in the anterior portion of 
the brain, had higher betweenness centrality for temporal - spatial 
(allpcorrected < 0.05, MFG: (1) = 9.90, IFTr: (1) = 35.90, IFOp: 
X^ (1) = 20.43, Fig. 4E). However, the MTL nodes showed levels of 
connectivity significantly greater than zero for both spatial and tem- 



poral order retrieval, indicating significant levels of connectivity dur- 
ing both forms of retrieval. Importantly, though, they exhibited no 
significant differences in their levels of connectivity for spatial versus 
temporal (all Pconected > 0.05, aHPC: x^ (1) = 4.30, pHPC: x' (1) = 
4.14, pPHG: x^ (1) = 0.11). Additionally, the node hierarchies in the 
right panels of Figure 4A,B emphasize a similar level of connectivity 
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Figure 4 | Differential Spatial and Temporal Network Connectivity. Panels A, B, and C are arranged similarly to Fig. 2 but show the spatial - temporal 
and temporal - spatial retrieval networks. (A, B) Connectivity plots, using standardized brain space (left plot) or hierarchical arrangement based on 
node degree (right plot), for spatial and temporal networks. (C) Left: Adjacency matrix where right of diagonal are significant nodes during spatial 
retrieval and left of diagonal are significant nodes during temporal order retrieval. Middle: Bar graphs comparing the number of edges for anterior versus 
posterior brain regions for the two conditions. Anterior brain regions showed higher connectivity during temporal retrieval while posterior brain regions 
showed higher connectivity during spatial retrieval. Right: Bar graph showing no difference in the total number of edges present in each network. 
(D and E) The MTL (aHPC, pHPC, pPHG) regions displayed high betweenness centrality values but no differences between condition. The SPL and 
calcarine ROIs showed significantly greater betweenness centrality for the spatial - temporal network, while IFG and MFG showed greater betweenness 
centrality in the temporal - spatial network. *Pcorrected < 0.05. 



and betweenness centrality for the different conditions. Contrary to 
our successful memory retrieval analysis, spatial and temporal 
retrieval networks exhibited similar density; a chi-squared test con- 
firmed that the number of connections is balanced between the spa- 



tial and the temporal conditions (x^ (1) = 1.38 p = 0.24, Fig. 4C, 
right), arguing that differences in regional connectivity could not be 
accounted for by lower degrees of connectivity in one condition. 
Again, differences in regional activation did not account for the 
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pattern of node connectivity as the two were uncorrelated (spatial - 
temporal: p = —0.17, p = 0.4120; temporal - spatial: p = —0.09, p = 
0.6740). Overall, these findings support the idea that spatial and 
temporal retrieval contribute somewhat independently to the correct 
retrieval effect we reported above, acting as subnetworks that overlap 
within the larger network involved in episodic memory retrieval. 

Discussion 

The present study aimed to test several hypotheses regarding the 
patterns of functional connectivity underlying retrieval of spatiotem- 
poral contextual details. By leveraging a new approach, whole brain 
functional connectivity with graph theory, we were able to better 
understand the contributions of different brain regions to episodic 
memory retrieval. Under the different conditions we examined (such 
as correct versus incorrect or spatial versus temporal trials), particu- 
lar networks emerged as important to episodic memory retrieval. 
Specifically, successful retrieval resulted in greater overall connectiv- 
ity combined with greater connectivity to specific hubs, including the 
hippocampus. Past fMRI studies employing univariate analyses have 
also focused on a primary role for the hippocampus in episodic 
memory^''''" although studies have also shown activation in the 
parahippocampal cortex'"^', prefrontal cortex"'^" precuneus"'"", 
visual cortex", posterior parietal cortex/ angular gyrus"''^, and retro- 
splenial cortex''"". These findings have not been unified across stud- 
ies, however, and univariate methods alone cannot indicate whether 
and to what extent these disparate areas interact"*. Thus, our results 
provide a potentially new perspective on the neural basis of episodic 
memory, emphasizing the importance of interactions across multiple 
brain regions rather than the computations of a single primary brain 
region (such as the hippocampus). 

Our findings comparing correct versus incorrect retrieval are con- 
sistent with several theoretical models of episodic memory"" which 
emphasize the importance of hippocampal-neocortical interactions. 
At the same time, our findings provide novel evidence that the hip- 
pocampus, and MTL more generally, indeed serves as a central "hub" 
for all forms of memory representation compared to other cortical 
regions. This is congruous with our previous findings using this same 
paradigm showing significant levels of activation in the hippocam- 
pus during both spatial and temporal order retrieval' and other 
studies emphasizing the importance of the hippocampus to both 
spatial and temporal order retrievaP \ The presence of the hip- 
pocampus as a primary hub in our study also provides novel insight 
into why the hippocampus frequently appears as a primary brain 
area in episodic memory. 

Our findings also bolster an important perspective in the memory 
literature: other brain hubs, outside of the MTL, contribute signifi- 
cantly to episodic memory. Past lesion work supports the idea that 
brain damage outside of the hippocampus, including prefrontal cor- 
tex'"" and parietal cortex'", impairs episodic memory, reinforcing the 
idea that this important function involves brain hubs across multiple 
lobes. Thus, whOe past work has also suggested the importance of 
extra-hippocampal brain areas to memory (for reviews, see"*''), our 
results suggest a reason why: prefrontal cortex, precuneus, and visual 
cortex are densely connected during successful memory retrieval. 
Our results further suggest one reason why hippocampal damage, 
in particular, may be devastating to episodic memory. While pre- 
frontal cortex showed greater connectivity during temporal order 
compared to spatial layout retrieval and parietal cortex showed 
greater connectivity during spatial layout retrieval compared to tem- 
poral order retrieval, only the hippocampus (and other parts of the 
MTL) acted as a hub during both spatial and temporal order proces- 
sing. The importance of the hippocampus to spatiotemporal memory 
is consistent with numerous studies using human neuroimaging''''\ 
intracranial EEG recordings'*, human lesion''*, rat electrophysio- 
logical studies'**, and rat lesion studies''"' suggesting its role in both 
spatial and temporal components of episodic memory. Thus, one 



possible interpretation of our results demonstrating important roles 
for the prefrontal cortex, precuneus, and visual cortex in episodic 
memory retrieval is that these areas likely facilitate and augment 
memory trace retrieval*''''^. One proposal, consistent with past work, 
is that prefrontal cortex and parietal cortex may be playing more 
auxiliary roles in memory, such as in executive or attentional pro- 
cesses'^'*^, with visual cortex supporting the visually-rich details that 
often accompany memory retrieval***. These more auxiliary roles are 
supported by the comparatively lower levels of connectivity they 
showed in contrast to the hippocampus in our study. By employing 
graph theory and functional connectivity in combination, future 
studies explicitly manipulating these variables, or meta-analyses of 
past studies, will be better able to address these possible specific roles. 

Our results converge with previous findings from our lab by 
Watrous and coUeagues""* who used a similar behavioral paradigm 
with intracranial EEG in human patients to show increased global 
connectivity across the frontal, medial temporal, and parietal lobes, 
with the MTL acting as a hub, during successful spatiotemporal 
retrieval. Previous findings have suggested similarities between 
low- frequency band oscillatory coherence and functional connectiv- 
ity measured using fMRI*', and corresponding results from our two 
studies (global increases in functional interactions for correct versus 
incorrect retrieval) are also consistent with the idea that low- 
frequency oscillatory coherence and function connectivity, measured 
via BOLD time series correlations, may carry similar types of 
information. Complementary to the current results, Watrous and 
colleagues found that spatial and temporal retrieval arose from com- 
mon anatomical but frequency-specific subnetworks, with somewhat 
limited regional specificity distinguishing between spatial and tem- 
poral order retrieval. While iEEG affords high spectrotemporal reso- 
lution compared to fMRI, fMRI provides potentially greater 
anatomical precision with greater ease of sampling a large number 
of brain hubs, including the hippocampus, which we did not sample 
in the Watrous et al. study. The two studies therefore highlight com- 
plementary aspects of retrieval, the spectrotemporal and functional/ 
anatomical substrates of memory retrieval. Importantly, our current 
findings thus provide several new important advances, implicating 1) 
the hippocampus as a hub in all aspects of spatiotemporal memory 
processing and retrieval 2) more comprehensive anatomical differ- 
ences for spatial versus temporal retrieval, particularly anterior ver- 
sus posterior differences, than previously observed 3) unique roles 
for the IFOr in incorrect versus correct retrieval and 4) greater levels 
of functional connectivity in individual participants as important to 
better memory performance. 

In conclusion, by employing a novel combination of fMRI, func- 
tional connectivity, and a graph theory approach, our study provides 
a new perspective on how the human brain processes episodic mem- 
ories. Our data provide novel support for models that emphasize 
global network interactions, rather than regionally mediated activity 
alone, as central to how we recover spatial and temporal memories 
associated with recent experiences. Our results nonetheless emphas- 
ize the importance of the hippocampus in particular as a convergence 
hub in this process, consistent with decades of research demonstrat- 
ing its importance to episodic memory. Our results thus suggest that 
while prefrontal and parietal areas communicate with the hippocam- 
pus during spatiotemporal memory retrieval, the hippocampus acts 
as a critical convergence "hub" during successful context retrieval. 

Methods 

Participants. Sixteen (13 female) healthy, right-handed participants participated in 
this study {mean age: 21.5, range: 20-28) and were compensated for their 
involvement. All were recruited from the University of California, Davis, and the 
surrounding communities and gave written informed consent, which was approved 
by the institutional review board at the University of California, Davis. Experiments 
and procedures were conducted in accordance with the IRB guidelines for 
experimental testing. 
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Experimental design and fMRI scanning procedure. During encoding, participants 
navigated a city consisting of eight stores randomly distributed around a circle in a 
virtual environment. They made a defined series of deliveries of a passenger to a store, 
with the spatial layout of the city uncorrected with the temporal order of deliveries to 
ensure dissociation between the contextual elements of an episode (Fig. lA). During 
retrieval, participants performed a spatiotemporal contextual retrieval task while 
being scanned. For each trial, participants were first asked whether they recognized a 
previously encountered store (item) compared to lure stores. To assay context 
memory, participants were presented with an additional question after each item 
recognition question regarding the contextual details associated with that item. Each 
trial was followed by a baseline task of variable duration based upon the jittered inter- 
stimulus interval (Fig. IB). The experimental design is further described in Ekstrom 
et al. (2011), where this dataset was previously analyzed using a univariate approach. 

The paradigm was configured as a mixed event -related, block design with 5 tem- 
poral and 5 spatial blocks that were interleaved for a total of 128 trials. Functional 
images were acquired on a Siemens 3 T TIM TRIO 32-channel scanner using a whole 
brain echo-planar imaging sequences (TR — 3 s, TE — 29 ms, slices — 47, field of 
view — 220 mm, flip angle — 90 degrees, bandwidth — 1684 Hz/pixel, voxel size — 2 
X 2 X 2 mm'). A high -resolution (voxel size — 1 X 1 X 1 mm') Tl-weighted 
magnetization- prepared rapid acquisition with gradient echo (MP-RAGE) structural 
image was also acquired to register functional activations to anatomical images. Using 
Statistical Parametric Mapping (SPM8) software (Wellcome Centre for 
Neuroimaging, UCL, London), the data were high-pass filtered, motion-corrected, 
and spatially smoothed with a 5 mm FWHM Gaussian kernel. Each individual par- 
ticipant's brain map was registered to the standard MNI template to carry out group- 
level analyses. 

Region of interest/anatomical parcellation. Rather than apply the traditional seed- 
based functional connectivity approach where seeds, or clusters of voxels, are 
determined a priori based upon known task involvement, regions of interest (ROIs) 
relevant to retrieval and other processes, such as attention and working memory, were 
sampled throughout the entire brain. Other functional connectivity methods, such as 
voxelwise comparisons of the entire brain'**, generate an unbiased, detailed profile of 
functional networks but encounter problems of multiple comparisons. The center of 
mass of twenty-six anatomical volumes delineated in standardized brain space'^" was 
used to determine the coordinates of thirteen bilateral ROIs. The ROIs selected for 
this study depended on the defined masks in the atlas, resulting in restricted 
availability for some regions, such as retrosplenial cortex. Each ROI consisted of a 10 
millimeter cube (a total of 125 voxels) assembled around the center of mass 
coordinate. If the constructed cubes encompassed voxels from surrounding regions, 
the cube was shifted in the x, y, or z coordinate directions to ensure the cube was 
comprised of voxels only belonging to the specified ROI. Finally, the hippocampus 
bilaterally was divided into anterior and posterior regions for a total of four ROIs, due 
to anatomical differences and functional specializations along the anterior to 
posterior axis^^ 

Statistical analysis and network construction. To assess inter-regional interactions 
throughout the brain, we employed a beta time series approach as described by 
Rissman and coUeagues^^. Briefly, each voxel's BOLD response in the task was 
modeled in a general linear model (GLM) as an individual regressor specifying the 
onset of each trial convolved with the canonical hemodynamic response function 
(HRF), resulting in 128 covariates of interest. The parameter, or beta, estimates 
derived for each trial for each voxel were then sorted by condition (spatial, temporal, 
correct, or incorrect trials) into a beta series. The spatial and temporal conditions 
consisted of both correct and incorrect retrieval questions to determine networks 
involved in the general process of spatiotemporal retrieval. The beta series of voxels 
belonging to an ROI were subsequently averaged culminating in twenty-six average 
beta series per condition. Finally, the Pearson product-moment correlation 
coefficient between all ROI's beta series was computed, creating a correlation matrix 
of all pairwise combinations (26 X 26) describing the strength of the functional 
relationship between two regions. 

To deal with the false positives, we derived a distribution of correlation coefficient 
values by correlating randomly shuffled beta series 1,000 times and pooling across 
participants and all pairwise ROIs within a condition. All correlation values were 
transformed using Fisher's Z, or the arc hyperbolic tangent transform^^. A functional 
connection, or edge, in the matrix was considered significant if the observed cor- 
relation value (correlation coefficient averaged across participants) was greater than 
the 90th percentile of the calculated distribution. This approach aimed to account for 
multiple comparisons and to fix the Type 1 error rate at 10%. This more liberal 
threshold preserves a higher network density thus providing enough connections to 
compute the graph theoretical metrics. Additional analyses using a bootstrap alpha 
level of 0.05 produced similar findings (Fig. S4) and suggest overall robust results. 
Additionally, the distribution of correlation values across all subjects and all pairwise 
regions was negatively skewed, so the average correlation value for all pairwise regions 
was greater than zero. Thus, we observed few negative correlations between regions, 
suggesting minimal "anti- correlated networks" within our findings. It is reasonable to 
expect that these areas might show mostly positive correlations in the first place 
because they have previously been implicated in memory processes. 

Those connections that withstood our thresholding procedure were assigned a 
value of 1 in the adjacency matrix, a binary, undirected matrix of pairwise connections 
between all ROIs in each condition. To compare conditions, the set intersection of 
spatial and temporal (spatial + temporal) and correct and incorrect (correct + 



incorrect) yielded elements that were common to both pairs of conditions (Fig. SI, 
82). In contrast, the set differences between the spatial and temporal conditions 
(spatial - temporal, temporal - spatial) and between the correct and incorrect con- 
ditions (correct - incorrect, incorrect - correct) revealed functional connections 
unique to a condition (Fig. 2, 4). 

Networks were then constructed in two-dimensional space, preserving the ana- 
tomical distance between brain areas in the x and y coordinate directions. In order to 
facilitate visualization of the connectivity across the brain, these two-dimensional 
networks were overlaid on top of a transverse slice of a standardized brain. Lastly, 
node hierarchies were created for each network where nodes are plotted vertically in 
ascending node degree order; nodes on the same horizontal level have the same node 
degree. For all graphs, the colored circles indicate ROIs, while a line between two ROIs 
signifies a functional connection. 

Behavioral measures. For the group level analyses, each individual participant's 
correlation coefficients were averaged to obtain a group correlation for each ROI pair, 
and then the group correlation matrix was thresholded as outlined above. For our 
behavioral analysis, we derived each individual participants network by thresholding 
that participant's correlation matrix again using bootstrap resampling, which resulted 
in 16 individual networks. The density of each participant's correct retrieval network 
was then computed and subsequently regressed against their behavioral accuracy for 
all task trials. Next, using stepwise regression, a final multiple linear model was 
produced to account for the variance in the density data. This final model was based 
on an initial model consisting of a subset of nodes whose percent contribution of the 
overall network density individually showed a significant (p < 0.05) linear 
relationship with performance across participants. The algorithm then stepped 
through all 26 regressors (or nodes) to construct a final model that accounted for the 
most variance (diagramed in Fig. 3). 

Network measures. The network topologies were characterized using measures of 
node degree, network density, and betweenness centrality on binary, undirected 
networks^'. Each measure was calculated from the Brain Connectivity Matlab toolbox 
(http://www.brain-connectivity-toolbox.net) or custom written Matlab code 
(Mathworks, Natick, MA). In the present study, a node, n, is defined as an ROI within 
the brain that has a defined anatomical location and volume within standard brain 
space, and N is the set of all nodes. An edge is defined as a significant functional 
connection between two nodes, where ajj is the connection status between node / and 
node j e N). ~ 1 when there is a significant connection, and ajj — 0 when no 
connection is present. Node degree quantifies the total number of edges connected to 
a particular node: kj — ^ a^. The density or total percent connectivity is the 
fraction of total number of edges in the network to the total number of possible edges: 
d — . ^ ' ^ . Finally, betweenness centrality is the sum of all the shortest paths 
that pass through the node of interest weighted by the inverse of the total number of 

EPhji'-) 
/jjew — - — , where p!,j is the shortest 

hi=!.h^m=i Phj 

path or the fewest number of edges that can be traversed between two nodes. For 
simplicity, we averaged the betweenness centrality across hemispheres (Fig. 2D and 
Fig. 4D,E). 
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